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Abstract 

It was recently demonstrated that static, resistive, magnetohydrodynamic 
equilibria, in the presence of spatially-uniform electrical conductivity do not 
exist in a torus under a standard set of assumed symmetries and boundary 
conditions. The difficulty, which goes away in the "periodic straight cylinder 
approximation," is associated with the necessarily non-vanishing character of 
the curl of the Lorentz force, j X B. Here, we ask if there exists a spatial profile 
of electrical conductivity that permits the existence of zero-flow, axisymmetric 
resistive equilibria in a torus, and answer the question in the affirmative. 
However, the physical properties of the conductivity profile are unusual (the 
conductivity cannot be constant on a magnetic surface, for example) and 
whether such equilibria are to be considered physically possible remains an 
open question. 
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I. INTRODUCTION 



Discussions of the magnetic confinement of plasmas and their stability usually begin 
with the subject of magnetohydrodynamic (MHD) equilibria.&i In the present generation of 
devices, dominated by the tokamak concept, the magnetofluids usually carry currents and 
contain strong dc magnetic fields whose sources are in external coils. 

Since the earliest days of fusion research, the various MHD equilibria have been treated 
as ideal. They are assumed to have infinite electrical conductivity and no flow, so that 
Ohm's law, 

E + vxB=j/a, (1) 

is satisfied by taking all three terms equal to zero. Here, E is the electric field, v is the fluid 
velocity, j is the electrical current density, and a is the electrical conductivity, allowed to 
become infinite. In familiar "Alfvenic" dimensionless units, the curl of the magnetic field 
B is just j, and of course, the divergences of B and j are zero. In the limit of infinite 
conductivity, the sharp connection that must exist between E and j for large but finite 
conductivity is, of course, destroyed. There have been attempts to reintroduce the effects of 
finite conductivity (e.g., see Grad and Hogani), but these attempts are in our opinion not a 
satisfactory resolution; too many possibilities are raised to be able to deal with any one of 
them conclusively. 

With the infinite conductivity assumption, the only remaining MHD requirement for 
equilibrium is one of mechanical force balance, with the local Lorentz force jxB taken 
equal to Vp, the mechanical pressure gradient. Such ideal equilibria are plentiful! and are 
not constrained by the Ohm's law. The Grad-Shafranov formulationi'iSi then provides a 
framework in which axisymmetric, toroidal ideal equilibria can be calculated. 

Investigating the stability of such ideal equilibria can be a lengthy undertaking. This is 
particularly true when the program includes investigating the various "resistive instabilities" 
that the ideal equilibria may be thought heir to. Aside from the questions of consistency 
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raised by ignoring resistive terms when constructing the equilibria, but then including them 
in the linear stability analysis, the possibilities for uncovering new variants of non-ideal, 
normal modes growing about ideal steady states can seem almost limitless. 

The purpose of this paper is to re-open the subject of zero-flow MHD steady states by 
considering the implications of retaining Eq. ([I]) with finite a. Attention will be confined 
to axisymmetric steady states with zero flow (v = 0) in the simplest toroidal geometry. 
A previous note by Montgomery and Shan! reported a proof of a somewhat unanticipated 
result: no such equilibria exist for the case of spatially uniform a. Here, we inquire into 
what kinds of spatial dependences of the conductivity o will permit axisymmetric, toroidal 
steady states without flow. We do not answer the question of how such conductivity profiles 
might occur, or indeed whether they are to be expected on physical grounds. 

The hydrodynamic precedent suggests the need for such inquiries, as illustrated by plane 
shear flows such as pipe flow, plane Poiseuille flow, or flow around obstacles (e.g., see 
AchesonB). In the case of pipe flow, for example, a uniform pressure gradient balances 
the axial viscous body force in the steady state. In the "ideal" limit, both may be allowed 
to go to zero proportionately, but the connection between them must not be lost, or what was 
a well-determined (parabolic) velocity profile may become anything at all; any axial velocity 
dependence upon radius is an acceptable steady state, clearly a nonsensical conclusion. Yet 
the stability of each one of these possible profiles could be investigated ad infinitum. 

It is perhaps worth remarking that in the absence of flow velocity, the situation de- 
scribed here is one of classical resistive magnetostatics, and may be treated by standard 
methods borrowed from electricity and magnetism textbooks. One can say that until a 
magnetofluid begins to flow, it does not know whether it is a solid piece of metal with a 
particular spatially-dependent electrical conductivity or a plasma. Once the requirement of 
balancing the Lorentz force against a scalar pressure gradient is met, none of the additional 
approximations or assumptions of magnetohydrodynamics per se need to be used for this 
problem. 

We also mention that recently, Bates and Lewis! have developed mathematical machinery 
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for treating toroidal equilibrium profiles in true toroidal coordinates. They were able to 
substantiate the result of Montgomery and Shani in considerably more detail. 

It is important that the point of view adopted here be stated clearly. We believe the 
full set of time-independent MHD equations, including resistivity, velocity fields, tensor 
transport coefficients, and realistic boundary conditions, are far too difficult to solve at 
this moment. This is the reason that nearly all MHD equilibrium theory has been ideal, 
and has omitted these effects, even when it has been asserted that "resistive instabilities" 
were being calculated. We are interested in moving as far as we can to remove the above 
restrictions, but we are also interested in doing so only by exhibiting explicit solutions to 
whatever problem is undertaken. We are fully aware of the desirability of including all of 
the complications. Here, only one is included: finite resistivity. It is a new result that there 
exist finite-resistivity toroidal MHD equilibria without flow, which, in many ways, do not 
look vastly different from ideal equilibria. How these states that are exhibited here relate to 
more general toroidal solutions with flow and the other possible complications is a separate, 
and much harder, question. 

This paper is organized as follows. In Sec. II, we ask and answer the question of what 
kinds of spatial profiles of scalar conductivity will permit axisymmetric toroidal equilibria in 
the presence of a stationary, toroidal inductive electric field. A partial differential equation 
is developed which can be thought of as replacing the Grad-Shafranov equation. In Sec. Ill, 
particular cases are worked out in detail. Implications are discussed in Sec. IV, along with 
questions about possible developments and extensions of the theory. 



II. THE FINITE-CONDUCTIVITY DIFFERENTIAL EQUATION 



In toroidal geometry (in contrast to the "straight cylinder" approximationtMliliy) ; the 
steady-state Maxwell equation, VxE = 0, has non-trivial implications for the zero-flow 
Ohm's law. We use cylindrical coordinates (r, <p, z) whose associated unit basis vectors are 
(e r , e^, e z ). The z axis is the axis of symmetry of the torus, and z = is the midplane. Ax- 
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isymmetry implies that the components of each vector field are independent of the azimuthal 
coordinate tp. The center line of the torus will be taken to be the circle r = r , z = 0. The 
toroidal cross section will be specified by giving some bounding curve in the r, z plane that 
encloses this center line. We will choose a specific shape when illustrating the results of the 
formalism. The wall of the container will be idealized as a rigid conductor coated with a 
thin layer of insulating dielectric. We do not consider the slits and slots that are necessary 
in any conducting boundary to permit the driving (toroidal) electric field to penetrate and 
drive the toroidal current. 

The electric field E = E^e^ will be assumed to result from a time-proportional magnetic 
flux (z direction) which is confined to an axisymmetric iron core whose axis coincides with the 
z-axis and which lies entirely inside the "hole in the doughnut" of the torus. The iron core is 
assumed to be very long, so the inductive electric field it produces in response to a magnetic 
flux that increases linearly with time is independent of the z-coordinate. Faraday's law then 
implies that E v (r) = E r /r only, where E is the value of the electric field at r = r . This 
is a highly simplified cartoon of a tokamak geometry, and it is not difficult to think of ways 
in which it could be made more realistic. For example, in a more refined model E might 
contain the gradient of a function of both r and z that obeys Laplace's equation. We do 
not consider here the possibility of a poloidal current density that could result from a tensor 
conductivity. 

The fact that E has this form, which is dictated by Maxwell's equations, has implications 
for j = erE: for a scalar conductivity a, j must also point in the tp direction. We shall assume 
that a is a function, as yet unspecified, of both r and z. 

The magnetic field consists of a toroidal part (B^ro/v)^, and a poloidal part B p (r, z) 
whose curl must be j: 



B is the strength of the toroidal field at r = r ; it is supported by external poloidal windings 
around the torus. The source for B p is j, so that 
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(2) 



VxB p =j = ^e,. (3) 

The boundary condition on the magnetic field is B ■ n = 0, where n is the unit normal to 
the surface of the torus. Thus the magnetic field lines will be assumed to be confined to the 
torus, and we do not inquire into the necessary external "vertical field," or other current 
distribution required to satisfy boundary conditions on B outside the torus. 

In addition to solving Eq. @ subject to the boundary condition, the remaining part of 
the problem is guaranteeing the force balance, Vp = jxB. This can be done by taking 
the divergence of this relation and solving the resulting Poisson equation for p. In order for 
there to exist such a solution, it is necessary that 

Vx(jxB) = Vx(j,e,xB p ) = 0. (4) 



For uniform a, it is not possibles to satisfy Eq. (f|). In the present circumstance, the question 
is asked: "What profiles of a(r, z) make it possible to satisfy Eq. @ and then ultimately, 
Eq. (0)?" Note that we are not necessarily assuming either incompressibility or uniform 
mass density. 

If we derive B p from a vector potential, B p = V xA^r, z)e v , and compute Vx(jxB), 
we see that Eq. (II) is equivalent to 



dA,d_ + l_d_ d_ 

dz dr r dr ^ dz 



^ = 0. (5) 



The general solution of Eq. (|||) is j v /r = firA^), where / is any differentiable function 
of its argument. The only remaining equation to satisfy is Eq. (||), or VxVxA^ = j^e^, 
which can be written as 

d 2 A ip 1 dA^ A v d 2 A ip 

-faT + r~dr + ~dz^ = ~ r/(r ^ ) - (6) 

Eq. (^|) is structurally similar to the Grad-Shafranov equation; any choice of / will lead 
to a partial differential equation for Atp(r,z) which must be solved subject to boundary 
conditions. However, a new restriction on the spatial conductivity is implied. Namely, the 
toroidal current density rfirA^) must determine a by 
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Eor a(r,2:) 

rf{rA v ) = • (7) 

The right hand side of Eq. fl7|) must always be positive, but seems otherwise unrestricted. 
In terms of the conventional poloidal flux function ip = rA^, the conductivity a is seen not 
to be a "flux function." That is, it is not constant on a surface of constant if), since it is the 
square of r times a function of ip. It is often thought that a should be constant on surfaces 
of constant if), since it depends strongly on temperature and flux surfaces are often imagined 
to be surfaces of constant temperature. 

For a specific form of the function /, the problem of finding the solution to Eq. (|6|) is 
straightforward; justifying the spatial dependence of the conductivity implied by Eq. (0), 
however, is not. 

We pass now to the consideration of some special choices of the function / and their 
consequences. 

III. TWO EXAMPLES 

We do not presently inquire into the physical basis of any possible choice of f(rA v ), such as, 
say, a maximum-entropy or minimum-energy choice. Rather, we select two examples largely 
on the basis of algebraic simplicity, and also choose a toroidal cross section for which the 
equation becomes separable and the boundary conditions tractable. Even then, both cases 
exhibit some complexity. 

A. The choice f = constant 

The simplest choice is to set / equal to a positive constant A > 0, making the toroidal 
current density j v proportional to r and independent of z. It should be noted that this 
r-proportionality (for vorticity) was found useful by NorburyEl in a study of vortex rings. 
These were ideal vortex rings, however, and required no external agency to maintain them 
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against dissipation. This makes the conductivity o vary as the square of r, increasing toward 
the "outboard" side of the torus. 

Eq. (P), with f — X, may be solved by first noting that it has become a linear, inhomoge- 
neous, partial differential equation, the general solution to which is any particular solution 
plus the most general solution of the associated homogeneous equation. A particular solu- 
tion is A v = — Ar 3 /8. The remaining homogeneous equation is the equation for a vacuum 
poloidal axisymmetric magnetic field. An easy way to find this vacuum field is to rewrite it 
in terms of a magnetic scalar potential $ that obeys Laplace's equation: 

B„ = VxA^ = -^Vxr\ + V$, (8) 

o 

where V 2 $ = 0. The boundary condition is B p ■ n = over the toroidal surface. 

Satisfying this boundary condition over a curved surface is a difficult task;! see, however, 
Bates and Lewis.! For illustrative purposes, we consider a torus with a rectangular cross 
section. We assume the toroidal boundaries to lie at z = ±L, and at r = r_ and r = r + , 
where r_ < r < r+. The vanishing of the normal component of B^ at these boundaries 
amounts to demanding that 

d$ . . 

— — = , at r = r_, r + (9a) 
or 

— = — , at z = ±L. (9b) 
The general solution of V 2( 3> = for $ can be written as 

$ = C z + £ [C k Jo(kr) + D k Y (kr)} sinh (kz + a k ) , (10) 

k 

where Jo and Y are the usual Bessel and Weber functions, respectively, and Co, Ck, Dj., 
and ttfc are arbitrary constants. The values of k remain to be determined. 
Eq. ( |9"aD is satisfied by requiring that 

C k J^kr_) + D k Y^(kr_) = 0, (11) 
C k J' {kr + ) + D k Y^(kr + ) = , (12) 



where the primes indicate differentiation with respect to the arguments of the functions. 
Eqs. (|Tl| ) and fll^D can only be solved consistently for C k and D k if the determinant 



V = J' (kr_)Y>(kr + ) - J',{kr + )Y^kr.) 

vanishes. For an infinite sequence of fc-values, with each k corresponding to a particular zero 
of T> for given values of r_ and r + , general Sturm-Liouville theory tells us that the functions 

<f) k = e k [J (kr) + (D fc /C fc ) Y (kr)} , (13) 

form a complete orthonormal set on the interval r_ < r < r + . The are real constants 
chosen to normalize the 4>h'- 



rdr = 5 kltk2 . (14) 



The ^-boundary conditions can both be satisfied by choosing the a k = for all k. The 
requirement is 

Ar 2 

C 4>o + ^2k C k 4> k {r) cosh kL = — , (15) 

where O — e o is a constant. This can be achieved if the expansion coefficients C k are chosen 
according to 

C = ^fJ r r+ r 3 dr, (16) 

for k — 0, and 



6* = ^ £ + fc (r)r 3 dr/ [k cosh fcL] , (17) 

for k 0. The full solution for B p is then given by Eq. with 

$ = Co 0o 2 + E Mr) smh(kz) . (18) 

In order to determine the allowed values of k in this problem, the zeros of the determinant 
T> must be found. The function D is an oscillating function of k that intersects the positive 
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fc-axis an infinite number of times. Using standard numerical techniques, the intersections 
can be computed for specified values of r_ and r + , and in this way the discrete spectrum of 
the permitted k values determined. For each k, we may calculate Ck/Dk. The results for the 
k values, the Ck/Dk ratios, and the normalization constants can be stored numerically, 
and the expansion coefficients Cj~ determined from Eqs. ( |T^ ) and (|TT[) . We may then plot 
magnetic surface contours, or surfaces of constant ip = rA^. The surfaces of constant 
pressure (since the current is purely toroidal) are guaranteed to coincide with the -^-surfaces. 

To see that contours of constant pressure coincide with constant ^-surfaces, it is useful 
to express the poloidal magnetic field B p in terms of ip: 

B p = V^xV^. (19) 

Then, the equation for scalar-pressure equilibrium, Vp = jxB, with j = Are^, can be 
integrated to give p = Xip + const. In Fig. 1, we illustrate surfaces of constant ip, for 
r_/r = 0.6, r + /r = 1.4 and L/r = 0.3. The current contours, in this case, will be rather 
strange, since the current is simply constant on lines of constant r, inside the torus, right 
up to the boundary. 



B. The case of linear variation 

The other case we wish to consider is that of a linear variation of /, or j proportional to the 
square of r times A v . The resulting linear differential equation is now homogeneous, which 
results in an interesting, but imperfectly understood "quantum" phenomenon: a preference 
for certain ratios of width to height for the rectangular cross section, in the steady states. 

It is convenient for this case to re-cast the differential equation in terms of the magnetic 
flux function if) = rA^: 

d fldip\ d 2 ip 

dr \r dr J dz 2 v 

We assume that the magnetofluid fills a torus with a rectangular cross section bounded by 

the planes z = ±L, and the radii r_ = ro — 5 and r + = ro + 5. Because the wall of the torus 
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is assumed to be perfectly conducting, we impose the boundary condition that the normal 

component of the magnetic field vanish at the wall. 

For j v = Xr ip, the equation for the flux function becomes 

d fldip\ d 2 ^ x 2 , n 
or \r or J oz A 

This equation can be solved in terms of confluent hypergeometric functions. To see this, we 
proceed as follows. First, we make the variable substitution p = \f~Xr 2 /2, which transforms 
Eq. (M) into 



Then, seeking separable solutions of the form ip(p, z) = R(p)Z(z), we find 



dz 2 

d 2 R / k 2 



- k 2 Z = 0, (21) 



v + [ 1 + «K) R=0 - (22) 

The equation for Z results in trigonometric or hyperbolic functions, depending on whether 
k is imaginary or real. Since the condition B p • n = on the rectangular wall of the torus 
leads to homogeneous boundary conditions on ip, 

^ = 0, at z = ±L (23a) 
Or 

dip 

—— = 0, at r = r + ,r_ (23b) 
oz 

the parameter k must be imaginary. Thus, the solution for Z is either sin^z) or cos(kz), 
where k = in, and k is real. The boundary condition in Eq. fl23a| ) can be fulfilled if we 
choose Z = cos(kz), and require k = n(n + \)/L, where n is an integer. 
The solution of Eq. (p2|) is given by 

R = C K Re [pe~ ip M(l - vq, 2; 2ip)] + D K Re [pe~ ip £/(l - ir], 2; 2ip)] , (24) 

where C K and D K are real constants, and rj = k 2 /A^/\. M(a,c;x) and U(a,c;x) are the 
regular and irregular, confluent hypergeometric (Kummer) functions, respectively^ They 
satisfy the second-order, ordinary differential equation 
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To fulfill the second boundary condition given in Eq. ( |23b| ), we demand 



C K Re 
C K Re 



e~ ip -M(l - ir], 2; 2ip_)\ + D K Re [e- ip -U(l - ir], 2; 2ipJ) 
e - ip +M{l - ir], 2; 2ip + )} + D K Re \e- ip+ U(l - ir], 2; 2ip + ) 



0. 
0. 



(25) 
(26) 



where p± = V\(r ± S) 2 /2. Eqs. (|25|) and (|26"D only have a solution for C K and _D K if the 
determinant vanishes, i.e., 



Re e~ lp - M{\ - ir],2;2ip. 
-Re \e- ip+ M(l - ir], 2; 2ip + ) 



x Re [e- ip+ U{l-ir],2;2ip^ 
x Re \e- ip -U(l-ir],2;2ipJ) 



0. 



(27) 



For each integer n, this equation holds only for a limited combination of A, r$, 5, and L. That 
is, for a particular value of the current density, only certain aspect ratios of the rectangular 
toroidal wall are allowed if steady state solutions are to exist. 

There is an additional physical constraint to be considered here which limits the values 
that the parameters n, A, r , 5, and L can assume: the current density must not change 
sign. Consequently, it seems that the only permissible value for n is n = 0. Note that 
the homogeneity of Eq. ( p0|) has eliminated the need for an eigenfunction expansion and 
satisfaction of boundary conditions, as in Eqs. ([16]) and (|T7|) . We have not experimented 
with the possibilities of linear combinations of the two types of / used in this section, but 
the range of possibilities is clearly wide. 

In the case that / varies linearly with its argument, the pressure is given by p = \ ip 2 /2 + 
const., which is easily obtained from integrating Vp = j X B, with j = Ar ip and B p written 
in the form of Eq. (p!9|). Thus, once again, the pressure does not vary on surfaces of constant 
-0 (magnetic surfaces). 

For n = 0, one set of parameters that satisfies Eq. ( p7| ) is Atq = 29.375, 5/r = 0.4, and 
L/r = 0.4. Contour plots of ip, j v , and a using these values appear in Fig. 2. One can 
see, from Fig. 2, that no radical qualitative departures from the magnetic surfaces expected 
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from a Grad-Shafranov treatment have been found; the "Shafranov shift" of the magnetic 
axis to the outboard side of the torus is clearly evident. The principal difference is the 
non-coincidence of the surfaces of constant conductivity and the magnetic surfaces. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have provided a framework in which toroidal axisymmetric, resistive 
steady states can been constructed for the MHD equations with scalar conductivity. Our 
approach is to search for a conductivity profile that permits such steady states to exist. 
Admittedly, the result is artificial since no discussion has been given of how such a profile 
might arise, or be consistent with an energy equation. If ideal MHD had not dominated 
the subject of magnetic confinement theory for forty years, the exercise might be considered 
unmotivated. However, since the proof of the nonexistence of toroidal resistive steady states 
with uniform conduct ivityi'i the question has arisen whether any current profile will support 
a static MHD state; here, we have answered this question affirmatively We consider the 
formalism for constructing static toroidal resistive states, incomplete as it is, to be physically 
less objectionable than an ideal treatment, which in fact still underlies the vast majority of 
instability calculations. 

What actually seems more likely, based on earlier dynamical computationJl^lliJii using 
the full three-dimensional MHD equations, is that velocity fields permit Ohm's law to be 
satisfied in resistive steady states of confined magnetofluids. Even if this conjecture is 
correct, though, the effects of finite flow will be bound up with that of the conductivity, 
which is virtually guaranteed to be spatially non-uniform. The satisfaction of the poloidal 
and toroidal components of both Ohm's law and the equation of motion, simultaneously, is 
an arduous task when velocity fields are included; at present, no one seems close to solving 
this problem. The situation becomes even more difficult if one demands that the pressure 
be derived from a local equation of state. In this case, complete consistency demands the 
simultaneous satisfaction of an energy equation as well - a formidably difficult undertaking 
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(e.g., see GoodmanEj). 

We should remark that the equilibria exhibited here could have been obtained formally 
as Grad-Shafranov equilibria with a constant poloidal current function, and with a proper 
choice of the pressure function. A demand for consistency with Ohm's law leads again to a 
determination of the conductivity through Eq. (0), with rA^ replaced by if), and / replaced 
by the derivative of the pressure function with respect to its argument (?/>). 

We should remark on a perception of several decades ago, due to Pfirsch and Schluter,0 
that difficulties associated with retaining Ohm's law and finite resistivity led to difficulties 
for ideal MHD equilibrium theory. Their resolution, though never absorbed to any significant 
degree in working models for tokamak equilibria, was to attempt to satisfy the Ohm's law in 
perturbation theory by using it to calculate iteratively the two perpendicular components of 
a fluid velocity to be associated with any given ideal equilibrium. A velocity field was thus 
taken into account in the Ohm's law, but not in the equation of motion. A perturbation 
expansion in inverse aspect ratio was also featured. The conclusion was reached that a 
universal plasma outflow existed (sometimes loosely called a "diffusion velocity") that had 
to be compensated by "sources" of new plasma within the plasma interior (never identified 
quantitatively). We cannot accept this conclusion. The explicit examples shown here, 
which do not rely on any perturbation expansions or other approximations, demonstrate 
the existence of zero-flow non-ideal solutions without plasma losses, and as such explicitly 
contradict the Pfirsch-Schluter expressions for positive-definite outward flux. We do believe, 
however, that Pfirsch and Schluter were correct in their assumption that real-life toroidal 
MHD steady states will involve mass flows in a fundamental way. It is in a sense remarkable 
that there has been so little attention to this as yet unresolved problem. 

The main point seems to us to be the need for developing a renewed respect for the 
problem of determining allowed steady states of a plasma. Neither present diagnostics nor 
present theoretical machinery permit it. The phasing out of the vocabulary that has arisen 
in connection with ideal steady states will require the passage of some time. 
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FIGURES 
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FIG. 1. Contours of constant ip in the case that j v = Xr, for r_/ro = 0.6, r + /ro = 1.4 and 
L/vq = 0.3. The ratio of magnetofluid pressure to poloidal magnetic field pressure (f3 p ) in this case 
is about 1.1. 
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FIG. 2. (a) Contours of constant tp (and pressure) in the case that j v = Xrip, for n = 0, 
Xtq = 29.375, 5/ro = 0.4, and L/tq = 0.4. (b) Current density contours, and (c) conductivity 
contours for the same set of parameters. The square box shows the location of the perfectly 
conducting, toroidal boundary. Note the non-coincidence of the conductivity contours and the 
magnetic surfaces. In this example, (3 P ~ 1.2. 
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